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Abstract. The cosmological evolution of an interacting scalar field model in which the scalar field interacts 
with dark matter, radiation, and baryon via Lorentz violation is investigated. We propose a model of 
£f>^ , interaction through the effective coupling /?. Using dynamical system analysis, we study the linear dynamics 

■ of an interacting model and show that the dynamics of critical points are completely controlled by two 

' parameters. Some results can be mentioned as follows. Firstly, the sequence of radiation, the dark matter, 

and the scalar field dark energy exist and baryons are sub dominant. Secondly, the model also allows 
the possibility of having a universe in the phantom phase with constant potential. Thirdly, the effective 
C^*) ' gravitational constant varies with respect to time through j3. In particular, we consider a simple case 

, where f3 has a quadratic form and has a good agreement with the modified /1CDM and quintessence 

q ■ models. Finally, we also calculate the first post-Newtonian parameters for our model. 

oo : 
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1 Introduction ory [6], the effect of a general class of such theories on the 

$-4 ' solar system has been considered in Ref. [7]. Moreover, 

d ■ There has been a growing appreciation of the importance as has been shown by authors in Ref . the Einstein- 

of the violations of Lorentz mvariance m recent years. The Aether theory may lead to signincant modifications of the 

intriguing possibility of the Lorentz violation is that an un- pQwer spectrum of tensor perturbation. The strong grav- 

known physics at high-energy scales could lead to a spon- itational cases including black holes of such theories have 

taneous breaking of Lorentz invariance by giving an ex- Deen studied in Refs [9] 
pectation value to certain non Standard Model fields that 

carry Lorentz indices, such as vectors, tensors, and gra- The existence of vector fields m a scalar- vector-tensor 

dients of scalar fields pQ. Recently, it has been proposed theor y of S ravlt y also leads to lts applications m mod- 

a relativistic theory of gravity where gravity is mediated ern cosmology and it might explain inflationary scenarios 

by a tensor, a vector, and a scalar field, thus called TeVeS UMM and accelerated expansion of the universe 0H2]- 

gravitational theory . It provides modified Newtonian The accelerated expansion and crossing of the phantom di- 

dynamics (MOND) and Newtonian limits in the weak field vlded lme has been studied recently by authors m Ref. [13] 

nonrelativistic limit, and is devoid of a causal propagation Based on a dynamical vector field coupled to the gravita- 

of perturbations. TeVeS could also explain the large-scale tl0n and scalar fields > we have studied to some extent the 

structure formation of the Universe without recurring to cosmological implications of a scalar- vector-tensor theory 

cold dark matter [3] , which is composed of very massive °^ g rav hy P3] • 

slowly moving and weakly interacting particles. On the Since the discovery of accelerated expansion of our 

other hand, the Einstein-Aether theory [4j is a vector- Universe |15j . identifying the contents of dark energy 

tensor theory, and TeVeS can be written as a vector-tensor and dark matter is one of the most important subjects 

theory which is the extension of the Einstein-Aether the- in modern cosmology. The dark energy is described by an 

ory [5]. In the case of generalized Einstein- Aether the- equation of state parameter u) — p/p, the ratio of the spa- 
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tially homogeneous dark energy pressure p to its energy 
density p. A value of u) < —1/3 is required for accelerated 
expansion. The classification of dark energy might be due 
to: quintessence field [16], tachyon models [17], Chaplygin 
gas rf5] if u) > —1, cosmological constant if w = — 1 [T9l 
|20"I|211I22] . or phantom field if oj < -1 [23]. A recent com- 
prehensive review on dark energy is available in [24] . Of 
course, as it has been discussed in [25 , 26 the vector field is 
also a viable dark energy candidate and effects on the cos- 
mic microwave background radiation and the large scale 
structure [2"T] . 

In the previous work [28] . the attractor solutions in 
Lorentz violating scalar-vector-tensor theory of gravity 
without interaction with background matter was stud- 
ied. In this framework, both the effective coupling and 
potential functions determine the stabilities of the fixed 
points. In the model, we considered the constants of slope 
of the effective coupling and potential functions which 
lead to the quadratic effective coupling with the (inverse) 
power-law potential. Differing from the previous work, 
in this work, we investigate the cosmological evolution 
of the scalar field dark energy and background perfect 
fluid by means of dynamical system. We study the cases 
of scalar field dark energy interacting with background 
perfect fluid. The interaction terms are taken to be two 
different forms which are mediated by the slope of the ef- 
fective coupling. For more realistic model we assume that 
the background matter fields might be dark matter, radi- 
ation, and baryons. 

Furthermore, to test the model in the solar system 
we present the post-Newtonian parameters (PN). In the 
PN approximation we restrict ourselves to the first post- 
Newtonian. The parameterized post-Newtonian (PPN) pa- 
rameters are determined by expanding the modified field 
equations in the metric perturbation. Then, we compare 
the solution to the PPN formalism in first PN approxima- 
tion proposed by Will and Nordtvedt [2"5]l30] and read off 
the coefficients (the PPN parameters) of post Newtonian 
potentials of the theory. 

This paper is organized as follows. In Section [2] we set 
down the general formalism of the scalar field interacting 
with background perfect fluid in the scalar-vector-tensor 
theory where the Lorentz symmetry is spontaneously bro- 
ken due to the unit-norm vector field. We derive the gov- 
erning equations of motion for the canonical Lagrangian 
of the scalar field. In Section [3l we study the interaction 
models and the attractor solutions. The critical points of 
the system and their stability are presented. The cosmo- 
logical implication is discussed in Section [4] In Section [5] 
we present the parameterized post-Newtonian parameters 
of the model. The final Section is devoted to the conclu- 
sions. 

In what follows, the conventions that we use through- 
out this work are the following: Greek letters represent 
spacetime indices, while Latin letters stand for spatial in- 
dices and repeated indices mean Einstein's summation. 
The symbol 0(N) stands for terms of order N. Finally, 
we use the metric signature (—,+,+,+). 



2 The action and field equations 

In the present section, we develop the general reconstruc- 
tion scheme for the scalar-vector-tensor gravitational the- 
ory. We will consider the properties of general four-dimensional 
universe, i.e. the universe where the four-dimensional space- 
time is allowed to contain any non-gravitational degree 
of freedom in the framework of Lorentz violating scalar- 
tensor-vector theory of gravity. Let us assume that the 
Lorentz symmetry is spontaneously broken by imposing 
the expectation values of a vector field as < 0| 1^1^10 >= 
— 1. The action can be written as the sum of four distinct 
parts: 



S — Sq + S u + S& + Sn 



(1) 



where the actions for the tensor field S g , the vector field 
S u , the scalar field S$, and the ordinary matter S m , re- 
spectively, are given by 



S a — 



S u 




Q 



ftVVV^ +A(u"u M + l) 
1 



(fx, 



-{V4>Y -v{4>) 



(2) 

(3) 
(4) 

(5) 



In the above f3i(4>) (i = 1,2,3) are the functions of 4>, 
and A is a Lagrange multiplier. In Eq. ([5]), we allow for 
an arbitrary coupling between the matter fields \P and the 
scalar field 4>. 

Varying the action ([1]) with respect to g^ v , we have 
field equations 



(6) 



where is the Ricci tensor, R is the scalar curvature, 
<? M „ is the metric tensor, and is the energy-momentum 

tensor for all the fields present, T M „ = T$ +T$ +T$ ) . 

Tfiy and are the energy-momentum tensors 

of vector, scalar fields, and ordinary matter, respectively, 
given by 

T$ = 2/3! (V„U T VuUr - V T u^ T u v ) - 2V r {u [p J T u) ) 
-2V T (u T J (/J „)) + 2V T (u (/1 Ju) T ) 
-2u a V T J T(7 u l _ l u u + g^Cu i (7) 

T$ = V^V„0 - \g^ [(V» 2 + 2V(4>)\ , (8) 



= (Pm +Pm)v fll V u +p m g fl u , 



(9) 



where is the four velocity and the current tensor 
in Eq. is given by 



(10) 
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The energy-momentum tensor is conserved 

v"(T#> + T$ + :r#>)=o. (n) 

In general, however, the Bianchi identity implies that each 
energy species in the cosmic mixture is not conserved, 
namely 

VT^ = , VT(« = o-W , V"rH = ff W (12) 

Here <r^ (fc = u, 0, m) is an arbitrary vector function 
of the space-time coordinates that determines the rate of 

transfer of energy, where cr^"' + a)jt + 4™' = 0. This is 
in accordance with Eq. (fTTj) . Equation (fl"2|) is the basic 
feature of interacting models in which there is exchange 
of energy between the components of the cosmic fluid. 
Moreover, the projection of the non conservation equation 
along the velocity of the whole fluid is 



Q 



-Q w -Q { 



(13) 



where = n^ajj,' is a scalar. 

Using Eq. ([6]), the Friedmann and Raychaudhuri equa- 
tions can be written as 



3H 2 = 8nG (p u + P4, + p m ) 



(14) 



and 



2-ff = SttG (p u + P0 + p m + Pu +P4>+ Pm) , (15) 

where 

-3/3H 2 , Pu = -p u + 2(pH + pH), (16) 



Pu 



1 



> 2 + V , P4 , = -P4, + 2 . (17) 



Here, we have defined /3 = (3\ + 3/32 + fh- 

Substituting Eqs. dTHJ) and (17|) into Eqs.CH]) and ( fl5 
respectively, we obtain 



(18) 



and 



2/3 + 



8ttG 



fl = -2/3H - 2 - 2(p TO + p m ) . (19) 



Let us define the effective coupling as follows 
then Eqs. (fT5)) and (fTi?|) can be simplified as 



H2 =<TP& 2 + V + P « 



H = _/3 Pm 

i? 2Hp lm H(3 



(20) 

(21) 
(22) 



Here, we have defined p m = (~f m — l) Pm , where 7 m is the 
ordinary matter barotropic parameter, which is related to 
the equation of state parameter oj m through the relation- 
ship 7 m = l+u> m . Similarly, we also defined the scalar field 
barotropic parameter, = (7^ — l)p<j, and 7^, = 1 + uj^. 
Then the effective equation of state for the total cosmic 
fluid is 



7 



(e) 



1 



Pu 



' Pn 



Pu + P<j> + Pn 



(23) 



which, again, is related to the equation of state parameter 
7^ e ) through 7^ = 1 + uj^ e \ The condition for an accel- 
erated universe is 7 (e) < 2/3. When < 7 (e) < 2/3, the 
universe is in quintessence phase while it is in phantom 
phase when 7^ < 0. 

From Eq. (|16[) we obtain 



(24) 



pu + 3H( Pu + Pu ) = 3H 2 P 



In order to preserve the conservation of total energy equa- 
tion p tot + 3H(ptot + Ptot) = 0, where p tot = p u + P4> + Pm 
and ptot — Pu + P<f> + Pm are the total energy density and 
the pressure, respectively, one can write the conservation 
of scalar field and matter field: 



P0 + 3H(p c/> 

in 



p<p) 

FPm 



= -3H 2 (5 + Q T 



(25) 
(26) 



The interaction term can be interpreted as a transfer from 
one energy component to another energy component of the 
cosmic fluid. These interactions are completely associated 
with Lorentz violation. In our case, the scalar field decays 
into the matter field and the vector field. The conserva- 
tion of scalar field, Eq. (j25|) . is equivalent to a dynamical 
equation for the scalar field cj>, 



Qr 



3H<p + V,* + 3H 2 p. 



(27) 



The above equation reduces to Refs. [Tl"ll2"5] for Q m = 0. 
Equations (f2"Tj) . (f2"2"| . and (|2T|) represent the basic set of 
equations of the model of interacting components of the 
cosmic fluid in the framework of Lorentz violating scalar- 
vector-tensor theory of gravity. In what follows we shall 
apply a dynamical system to analyze the cosmological dy- 
namics of this set of equations. 



3 Interacting model 

Some models that allow interaction between the scalar 
field and the matter field have been proposed as a solu- 
tion to the cosmic coincidence problem. These models are 
compatible with observational data but so far there has 
been no evidence on the existence of this interaction. A 
solution will be achieved if the dynamical system presents 
scaling solutions which are characterized by a constant 
dark matter to dark energy ratio. Even more important 
are those scaling solutions that are also attractors and 
have the accelerated solution. In this way, the coincidence 
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problem gets substantially alleviated because, regardless 
of the initial conditions, the system evolves towards a fi- 
nal state where the ratio of dark matter to dark energy 
remains constant. 

The explicit form of Eq. (fT3|) can be expressed in the 
form 



= -^{P<P+ Pm) 



We assume the interaction term as follows 



(28) 



n -P 

P 



(29) 



The interaction term (|29|) means that the scalar field can 
exchange energy with the background matter, through the 
interaction between them. In this case the exchange energy 
is mediated by the slope of the effective coupling. 
Equations (|2"5|) and (f2l)]) . respectively, become 



P4, + 3H(p <j> +P4>) 

H in 

+ 3H(p 

m t Pm 



P± 

P 

P. 



P 



(30) 
(31) 



For a more realistic model we assume that the matter 
fields might be a combination of dark matter, p c , radia- 
tion, p r , and baryons, pb'- p m = Pc+Pr+Pb- We also assume 
that the barotropic equation of state for the radiation field 
Pr = Pr/3 and that the baryons are non-relativistic par- 
ticles so that pb = holds. Hence, the equations for the 
energy densities of radiation and baryons are 

p r + AHp r = , p b + 3Hp b = , (32) 

respectively, and we find the well-known relationships: 
p r oc a~ 4 and pb oc a" 3 , a is a scale factor. For the scalar 
field and the dark matter we have 

p <j> + 3H^ ) p <j> = , p c + ZHjMpc = , (33) 

where 7^ and 7c are the effective barotropic equation 
of state for scalar field and dark matter, respectively, 



70 



3H/3 



P 



3H/3 



Notice that for P/P < we have 7^ < 70, r y c c> < j c and 
both and p c with Lorentz violation will dilute slower 
then that without Lorentz violation or /3 = const. Thus 

P/P will determine both the effective equations of state 

(<0 i (e) 

7^ and 7c . 



Pr + Pb 
Pc 



» 



(34) 



3.1 Dynamical analysis 

In order to study the dynamics of the model, we shall 
introduce the following dimensionless variables [T41I28] : 



2 _ 



6#ff 2 

P,0 



y 



v 



3H 2 p 



(35) 
(36) 



i l = ^5 — > l2 = — T79 1- 77 t^ nr ) (37) 



0. 



vl 2V,*/v 



and, accordingly, the governing equations of motion could 
be reexpressed as the following system of equations: 

H' = -^h(i + x 2 - V 2 + ^z 2 -V6X 1 x\ , (38) 
H 1 

■>■' = - < (3 + 77 



+ \Jl (Ai+A 2 )y 2 + 2^A 1 x 2 



where 



3/3iA 



(39) 
(40) 

(41) 

(42) 

(43) 



A prime denotes a derivative with respect to the natural 
logarithm of the scale factor, d/d\na = H~ 1 d/dt. Equa- 
tion (f2"Tj) gives the following constraint equation: 



fir 



Pc 



3/3H 2 



= 1 - x & 



(44) 



where = p$/3PH 2 
and Q b = p b /3f3H 2 = 



y 2 , Or = Pr/3~PH 2 



Notice that fli, (i = (f>,c,r,b) 



are the effective cosmological density parameters which 
are associated with the Lorentz violation. In general, the 
parameters Ai, A2, A and A are variables dependent on 
<fi and completely associated with the Lorentz violation. 
In order to construct viable Lorentz violation model, we 
require that the effective coupling (3 and the potential 
function V should satisfy the conditions A > 1/2 and 
r% > 1 — A1/2A2, respectively. In this paper, we want to 
discuss the phase space, then we need certain constraints 
on the effective coupling and potential function. Note that 
for Pi = const., Ai — > 0, the scalar field dynamics in the 
Lorentz violating scalar-vector-tensor theories is then re- 
duced to the scalar field dynamics in the conventional 
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Table 1. Properties of the critical points. 



Point (x,y,z,u) 



7<A 



7e// 



A- 
B 

C. r 

D 

D r 

E b 



(+1,0,0,0) 
(-1,0,0,0) 




,0,0 



2-4/2/3 / 4(A 2 -2Ai) 
*2 ' V 3A|(A!+A 2 ) ' 

(0,0,1,0) 
(0,0,0,0 



(4A 1 +A 2 )(A 1 A 2 + A|-4) 



1 
1 
1 

2 
3 A j 

3 

(Ai+A 2 )2 
4 

A 2 (A!+A 2 ) 






2 
2 

(Ai + A 2 ) 2 
3 



4(A!+A 2 ) 
3A 2 



2 - Ai 
2 + 2,/|Ai 



2 

3 

1 _ 2A 1 
Ai + A 2 

4 / X2-A1 



3 A 2 



Table 2. 

points. 



Stabilities and acceleration conditions of the critical equation of motion 



Point 


Existence 




Stability 


Acceleration 


A+ 


VAi, A2 




unstable 




A. 


VAi, A2 




unstable 


A i<-\/i 


B 


(Ai + A 2 ) 2 


< 6 


stable 


A| < A? + 2 


Cf 


Af >| 




unstable 


never 


D 


(Ai +'A 2 ) 2 


> 3 


stable 


A 2 < 5Ai 


D r 


A 2 (Ai + A 2 


)>4 


unstable 


A 2 < 2Ai 


Ey 


VAi, A2 




unstable 


never 


E b 


VAi, A 2 




unstable 


never 



one. But, the effective gravitational constant is rescaled 
by Eq. (J2DJ). In this case, the cosmological attractor solu- 
tions can be studied through a scalar exponential potential 
of the form V{(f>) = Vq exp(— A20/ where (3 = const. 
This exponential potential gives rise to scaling solutions 
for the scalar field [51] . In this paper we consider the case 
in which Ai and A2 are constant parameters. For example, 
a constant Ai is given by an effective coupling f3 — £0 2 and 
we have Ai = —2y^._A constant A2 can only be obtained 
as a combination of (3(4>) and V(<j>), one finds 



v{cj>) = v {M.4>)Y , 



(45) 



where s — A2/A1 is a constant parameter. In general, one 
can write the potential as a function of effective coupling, 
= /(fly)). 



3.2 Attractor solutions 

The critical points (x c , y c , z c , u c ) are obtained by imposing 
the conditions x' = y' = z' = v! = 0. Substituting linear 
perturbation x — > x c + 5x, y — > y c + Sy, z — > z c + 5z 
and u — * u c + Su about the critical points into Eqs. ([55]) - 
(|42p . we obtain, up to first-order in the perturbation, the 



\V\=M 





(46) 



Notice from ([39]) -([42 | that the dynamical equations are in- 
variant under the change of sign (y, z, u) — * (— y, —z, —u), 
and in consequence we dont have to include the points 
with (y, z, u) < in our analyzes. The properties of the 
critical points are summarized in Table [T] There are eight 
critical points at all and two of them lead to attractor 
solutions, depending on the values of the parameters Ai 
and A2. The scalar field dominated solution, point B in 
Table [5J are characterized by Q = 1 , and the effective 
equations of state are given by 

1 



r (Ai + A 2 ) 2 , 7 (e) 



1 



(A? 



A) 



(47) 



The solution of this point exists for (Ai + A2) 2 < 6 
and the universe is accelerated for A?, < \\ + 2. From 
eq. ([4"T]) one can see that the de Sitter epoch corresponds 
to A2 = Ai. The scalar field is dark energy when A 2 < 1/2. 
In this case the effective coupling (3 and the potential func- 
tion are quadratic in <fi, j3(<f>) ~ V((f>) ~ 4> 2 . The inflation- 
ary solution of this model has been studied in Ref. 
Figure [1] shows that the sequence of radiation, dark matter 
and scalar field dark energy. The baryon is sub-dominant 
in this case. The parameters correspond to A2 = Ai and 
Ai = — l/y/3. The scalar field equation of state parameter 
w — l<p ~ 1 is nearly a constant, during the radiation 
and matter epochs because the fields are almost frozen 
for which lo^ = oA e ^ . At the transition era from matter 
domination to the scalar field dark energy domination, 

(e) 

w<j, and uj)p begin to grow because the kinetic energies of 
the fields become important. However, the universe enters 
the de Sitter phase during which the field 4> rolls up the 
potential. More interesting of this attractor solution is of 
the constant potential, A2 = 0. The universe is in phan- 
tom phase in this case because uj^ is crossing —1 and it 
is accelerated for A 2 > —2. 
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Fig. 1. Evolution of the density parameters and the equation of state parameters as a function of lna. Top panel corresponds 
to the case of A2 = Ai = — l/v3 while the bottom panel corresponds to the cases of constant potential and Ai = — 1. 
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Fig. 2. Evolution of the density parameters and the equation of state parameters as a function of lna. Top panel corresponds 
to the case of A2 = Ai = — 3/y2 while the bottom panel corresponds to the cases of constant potential and Ai = 3/2\/2. 



The second attractor solution is the scalar field scaling 
solution, point D in Table [2] The solution of this point 
exists for (Ai + A2) 2 > 3, corresponding to energy density 
parameter = 3/(Ai + A2) 2 . The effective equations of 
state are given by 



14> = In 



„( e ) 



v( e ) 



A 2 



(A1+A2) 



2Ai 



Ai + A2 



(48) 
(49) 



The universe is accelerated for A2 < 5Ai. In the case of 
the effective coupling [3 and the potential function are 
quadratic in 0, i.e. A2 = Ai, the universe is always accel- 
erated. For the constant potential, A2 = 0, the scalar field 
behaves as a cosmological constant while the universe is in 
phantom phase. Figure [H shows the sequence of radiation, 
dark matter and scalar field dark energy. The baryon is 
sub dominant in this case. The parameters correspond to 
A2 = Ai = — 3/V2 (top panel), and Ai = 3/2\/2 (bottom 
panel) . 
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4 A comparison of the model using supernova 
data 

From the above detail analysis, we may investigate the 
cosmological consequences of a Lorentz violating scalar- 
vector-tensor theory which incorporates time variations 
in the gravitational constant. It was raised by Dirac who 
introduced the large number hypothesis [32], and has re- 
cently become a subject of intensive experimental and the- 
oretical studies [32]. The effective gravitational constant, 
G^ , is obtained from the Friedmann equation, 



G (e) = -L = G 

8tt/3 1 + 8ttG/3 ' 



(50) 



From Eq. (|55| . it is clear that the Lorentz violation leads 
to time variation of the gravitational constant. 

In the following, we study the expansion history of the 
universe using the 194 Snla data [35]l36], We simplify our 
model by considering an interaction between dark matter 
and the scalar field dark energy given by Eqs. (|30]) and 
([3Tjl . The evolution of the dark matter and scalar field 
dark energy are given by 

pt(z) = Pi0 e 3 ^ W , (i = m, 0) , (56) 

where z = 1/a— 1 is the redshift. Using the above relation, 
the Hubble parameter as a function of the redshift can be 
written as 



where G is the parameter in the action ([T]) . Therefore the 
time variation of G^ can be written as 



G(e) 







(51) 



and the effective gravitational constant is determined by 
the effective coupling /?. For the quadratic effective cou- 
pling, {3 oc 2 , the effective gravitational constant is in- 
versely proportional to 4> 2 , G^ oc [cf>(t)]~ 2 . Recently using 
the data provided by the pulsating white dwarf star G117- 
B15A the astereoseismological bound on G/G is found [34] 
to be -2.5 x 1CT 10 yr- 1 < G/G < 4.0 x 10" 10 yr' 1 . In 
the present model the time variation in the gravitational 
constant is given by 



G^ 



3Ai 



G( £ ) (Ai + A 2 ) 
in the scaling solution and 



H 



G( £ ) 



Ai(Ai + A 2 )ff 



(52) 



(53) 



in the scalar field dominated solution, where the evolution 
of the Hubble parameter is given by Eq. (|38|) . For instance, 
in the case of power law expansion of the universe a(t) oc t p 
with p > 0, the time variation of G^ leads to 



3Ai 



G( £ ) (Ai + A 2 ) 



(54) 



in the scaling solution. Assuming the present age of the 
Universe as 14 Gyr, it is straightforward to derive from 
Eq. (H]) the estimate ~ 2.14 x 10" 10 yr~ x for 

the case of constant potential. Our model also allows the 
negative value of G^'jG^ 1 . Let us focus on the scaling 
solution. If Q$ — 2/3 we find 



G (e) 
GW 



= ±V2XiH 



(55) 



A negative G^ / G^ implies a time-decreasing G^ , while 
a positive G<» /G^ means is growing with time. 



+ (l-r? m0 )(l + z) 3(1+ ^ (z)) ] , (57) 

where the subscript describes the current value of the 
variable. Notice that the evolution of the Hubble param- 
eter is deviated by the factor of_(/?o//?) 2 compared to the 
standard one. If the functions 0(z) and ui^z) are given, 
we can find the evolution of the Hubble parameter. In this 
section, we consider an ansatz for the effective coupling, 



£ = #,(! + O 2 ) 



(58) 



where £ is a constant. 

Let us first consider the modified A Cold Dark Matter 
(ACDM) model. We have 



H 



\z; C, n m0 ) = (7-^2) [ fl mo(l + z ? + (! - ^™ol59) 



Equation (|59p has two free parameters £ and f2 m o which 
are determined by minimizing 



^2 _ \pobs(zi) - fi(zi)Y 



a. 



(60) 



where \x is the extinction-corrected distance modulus, 



^ = 5l0g -(l^) +25 ' 



(61) 



and tji is the total uncertainty in the Snla data. The lu- 
minosity distance is given by 



d L (z) 



c(l + z) f z dz 



Ho Jo H{z>) " 



(62) 



Fitting the model to 194 Snla data, we get Xmin = 
195.68, C = —0.33, and fl m o = 0.24. For comparison, we 
also fit the cosmological constant model to the 194 Snla 
data and find \ 2 = 198.74, and Q m0 = 0.34. 

In the next model we replace the cosmological constant 
energy density by a scalar field dark energy with constant 
equation of state parameter (lo^,(z) — constant). We set 
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0.25 0.5 0.75 1 1.25 1.5 1.75 0.25 0.5 0.75 1 1.25 1.5 1.75 



z z 

Fig. 3. Observational 194 Snla Hubble free luminosity distances fitted to our model. Left panel corresponds to the case of the 
cosmological constant. The best fit values are £ = —0.33, fi m o = 0.24. Right panel corresponds to the case of the quintessence 
with constant equation of state parameter. The best fit values are C, = —0.29, ui^, = —1.13. Continuous line denotes the curve 
in the context of Lorentz violating scalar-vector-tensor theory, while dashed line denotes the standard one. 



here f2 m o — 0.3. We evaluate X 2 {C,i u3 <t') an d minimize with 
respect to ( and We find 



X 2 (C = -0.29, uj^ = -1.13) = 195.71 . (63) 



Figure [3] shows a comparison of the observed 194 
Snla Hubble free luminosity distances along the predicted 
curves in the context of Lorentz violating scalar-vector- 
tensor theory. We see that the effect of Lorentz violation 
appears at z > 0.75. We define the reduced form of Hubble 
parameter compared to the standard case as 



tj2 _ tt2 
tt2 _ n LV n std 
n red tj2 



(64) 



where 



n 2 

n std 



+(i-n m0 )(i + zf 1+ ^ z » 



(65) 



Thus the reduced form of Hubble parameter, due to the 
effect of Lorentz violation, is 



Hred( z ) 



0(z) 



5 Parameterized Post- Newtonian 

In order to confront the predictions of a given gravity the- 
ory with experiment in the solar system, it is necessary 
to compute its PPN parameters. The post-Newtonian ap- 
proximation is based on the assumptions of weak grav- 
itational fields and slow motions. It provides a way to 
estimate general relativistic effects in the fully nonlinear 
evolution stage of the large scale cosmic structures. The 
procedure of parameterizing our model is following to that 
of Ref [37] , in which the authors has derived the PPN pa- 
rameters in the frame of Einstein aether theory. 

In the weak field approximation, we choose a system of 
coordinates in which the metric can be perturbatively ex- 
panded around Minkowski spacetime. The decomposition 
is as follows 



where rj^ is the Minkowski metric and h^ v is the metric 
perturbations and we take |/i M „| << 1. 

The equations governing the perturbation h^ lL , in the 
model Eq. (p} are found by computing the Einstein field 
equations in the perturbative limit. The full field equations 
are given by 



We allow for an arbitrary coupling between the matter 
fields and the scalar field 0. We assume that the scalar field 
</> is coupled to a barotropic perfect fluid with a coupling 
function given by 



1 



5S„ 



(69) 



Pm\j-g Of 

Therefore, the equation of motion for the scalar field is 



(70) 



(66) where 



K x = VVV^, K 2 = (V/f , K 3 = Wu"V u u£71) 

In the previous discussion we have considered the coupling 
between the scalar field and the matter field is given by 
the effective coupling, q(4>) — dln0/d<p, where (3 is defined 
by Eq. (HU). 

We also have the vector field equation, 



(72) 



where J M „ is given by Eq. (fT0|) . The constraint for the 
vector field is 



g^u v + 1 = 



(73) 



The vector field is purely timelike at the zeroth order and 
the fluid variables are assigned orders oip^II^p/p^ 
(v 1 ) 2 ~ 0(1). The scalar field is expanded as 



(67) 0(1), where 4>q is determined by the cosmological solution. 
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Then, the metric perturbations h^ v will be of orders hoo ~ 
0(1) + 0(2), hj ~ 0(1), and h oi ~ 0(1.5). 

The general form of the first post Newtonian metric is 
given by [30 



9oo 



9oi 



-1 + 2U- 2(3 PPN U 2 - 2f#w - (Ci - 20-4 
+ (2 lPPN + 2 + a 3 + Ci - 20*1 
+2{i lPPN - 2/3 PPN + I + C2 + 0*2 
+2(1 + Cs)*3 + 2(3 lPPN + 3( 4 - 20*4 , 
(1 + 2 1 p PN U)8 ij , 

- g {^Ippn + 3 + ai - a 2 + Ci - 20K 



-2(1 + 02 -Ci- 20 Wi 



(74) 



The PPN potentials (J7, <? w , *i, *2, *3, *4, A Vi, Wi) are 
defined by 



F(ar) = G N / d J y 



I I s 

F - V\ 

where the correspondences F : f are given by 
U : 1, #1 : 2 : ?7, *3 : *4 



(75) 



P 



<f>„ : / d 3 z Pm (z)^U ( ^-4*Y r, : r'. 



.4 



k - y| 2 V k ~ A ) ' 
]r :Vj{x - y)j{x - y y 



\x-y\ 2 ' ' l^-yl 2 
These potentials satisfy the following relations 

F u = -ATrG NPm f , (77) 
for U, *i,2,3,4) and V,. The superpotential x is defined by 



X = -G N / d 3 yp m \x - y\ , 



(78) 



which satisfies 



X,» = -2U 
We also note the identity 



X,o* =Vi-Wi 



(80) 



lowest Newtonian order. One can verify that in general 
relativity jppn — Pppn — 1 an d all other parameters 
vanish, which implies that there are no preferred-location 
or frame effects and that the theory is fully conservative. 

The PPN parameters are determined as follows: ex- 
pand the modified field equations in the metric perturba- 
tion and in the PN approximation; iteratively solve for the 
metric perturbation to 0(2) in /loo, to 0(1.5) in and 
to 0(1) in hij] compare the solution to the PPN metric of 
Eq. (|74|) and read off the PPN parameters of the theory. 

The gauge condition we use for the metric is as fol- 
lows [37] 



hot 



3£/ + Bu\ , 



(81) 
(82) 



where B is a function of Pi which will be determined be- 
low. 

There are two important notes to be considered here 
since we concern in the post-Newtonian expansion. Firstly, 
the term V(4>o) is the same order as the energy density 
of the cosmological constant, therefore, these term can- 
not leads to any observable deviations at Solar system 
scales. Practically, we can assume dV (cf>) / d<fi — as far 
as the post-Newtonian expansion is concerned. It means 
that the cosmological solution corresponds to a minimum 
of the potential and the scalar field will satisfy the ex- 
ternum condition. Secondly, from the analysis of the ten- 
sor perturbations, the velocity of the gravitational waves 
are different from the velocity of light [11 . In order to 
have the real velocity, the coupling functions Pi and P3 
must be satisfy (Pi + Pz) < (167rG) _1 . Therefore, one can 
assume that pi and P3 are constant. This assumption does 
not affect our previous calculations because the coupling- 
functions Pi and pa have included in the definition of the 
effective coupling, which is given by Eq. ptJ|) . These two 
assumptions will be used in the following calculations. 



(79) 5.1 Solving u^, <f), and g^ to 0(1) 



We first solve the constraint equation ([73")) to 0(1) 
find 



The PPN metric Eq. ifTi]) contains ten parameters (jppn, 
Pppn, a i.2,3j Ci.2,3.4). The parameter 7ppjv measures 
how much space-curvature is produced by a unit rest mass, 
while the parameter Pppn determines how much non- 
linearity is there in the superposition law of gravity. On 
the other hand, the parameter £ determines whether there 
are preferred-location effects, while 01,2.3 represent preferred- 
frame effects. Finally, the parameters £1,2,3,4 measure the 
amount of violation of conservation of total momentum. In 
terms of conservation laws, one can interpret these param- 
eters as a measure whether a theory is fully conservative, 
i.e. the linear and angular momenta are conserved ((1,2.3,4 
and 01,2.3 vanish), semi-conservative, i.e. the linear mo- 
mentum is conserved ((1,2,3,4 and 0(1,2,3 vanish), or non- 
conservative, where only the energy is conserved through 



u u = 1 

From this result, we have 
1 



\hoo 



U = -1 



where we used u„, = 



'oh 



hoi 



From Eq. (|TU)) . we obtain 



qoPr, 



= . 



Using Uji — —AirGNPm, we have 



4nG 



-U A 



N 



We 



(83) 



(84) 



(85) 



(86) 
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The solution of this equation is 



« = S°-U = -^-U, (87) 



AttG 



N 



4ttG 



N 



where qo — q(4>o) and ln'/3o = dhi /3(<po) /d<j>. 

The general expression of the covariant derivatives of 



VX^^ + yrtV^ + W-o-V,") • (88) 
To 0(2), we find 

V M u° = 0, \7 u l = Uj, - ^h 00: i + ^h ohoo,i , (89) 

and 

11 3 

Vi(V u l ) = tiifii - -hoo,u + -^hoohoo.u - -^h 00ii hoo,i{90) 

We also have 

ViU J = Uj,i + - (hjifl - h 0j .i - h i,j) , (91) 
to 0(1.5). 

Let us now proceed with the solution to the field equa- 
tion for the "time-time" component of the metric pertur- 
bation. The left hand side of Eq. (|69]) is 



Roo — —-zhoo,u + -zhijhoo,ij + ( hiQi — —hu,o 



1 



to 0(2). Thus, at 0(1) we have 

■Roo = — g^oo.M • 



1 



(92) 



(93) 



To 0(1), we have the components of energy-momentum 
tensor 

T^= Pm , Tj l) =0, T$=0 = T™, (94) 

and 

Tqq' 1 = — 2VoJoo — 2V M Jo M = — 2VjJo l 

= 2(/3iV u l ) = -/3io V« - (95) 

(96) 

here fiio = A(<Ao), i = 1, 2, 3. Then, we obtain 

(1 - 8ttG/3io)Jioo,« = -8irG Pm , (97) 
which gives hoo to 0(1), 

/ioo = 2U , (98) 
with Newton's constant 

G 



rp(u) _ j-. 



G 



JV 



(1 - 87rG/3 10 ) ' 



(99) 



To 0(1) the "space-space" component of left hand side 
Eq. S is 



Rij — —-zhij^kk — -zhkksj + 7:^00, ij 



If we use the gauge (l8~Tj) . we obtain 



ft 



1 



r/l 



ij,kk 



(100) 



(101) 



Then, we have 

(1 - %TTGp X0 )h ijM = SwGpmSij . (102) 

This expression similar to the "time-time" component. It 
means that the spatial metric perturbation to 0(1) is then 
simply given by the GR prediction without correction, 
namely 



where hoo is given by Eq. ([98]). 



(103) 



5.2 Solving u\g 0l to 0(1.5) 

From the previous results one can see that the Lagrange 
multiplier is A ~ 0(1)- Therefore, the vector field equation 
yields V M J^i = 0. This means 



Joi,o + Jji.j — , 



(104) 



to 0(1.5), where 



Joifi = -C3iV Ui),o , (105) 
Jji,j = -(/3iV,-Ut)j - (^ 2 V fe u fe ) : i - (feVjU^j (106) 

We notice that V it l = -\h Q0 ,i to 0(1). To 0(1.5), we 
have 



^0i,0 — ^/^lO^OO, iO j 



(107) 



and 



■I,, j = -PioUjj - (P20 + P3o)u 3 yij ~ (010 ~ P30)h [i,j]j 
-\(fiia + aP2D + Pso)hao,oi- (108) 
Then, the vector field equation can be written 

PlOUjj + (#20 + P30)u% + (/3l0 - 03o)ho[i,j]j 

-|(2/?io + 3/3 20 + Pao)x,oiji = . (109) 

By taking the spatial divergence of this equation, we can 
solve for u % t . 



(110) 
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where 



C = 



(2/3io + W20 + &o) 



2(/3io+/3 2 o+/3 3 o) 
Using Eq. (|110|) , the gauge (|82p can be rewritten 
1 



hoi 



-(3 - 2£G )x, 0i 



(111) 



(112) 



Substituting Eqs. pTU|) and (TT1D into (HOHD a and us- 
ing the previous results, we can solve Eq. (|109|) for u l , 



(0w ~ Po) ■ 
2/3io 

Co 



(/3 10 -A))(3-2BCo) 
4/3io 



X,oi 



(113) 



Let us now look for solutions to the field equations for 
the metric perturbation g 0i to 0(1.5). To 0(1.5), we have 
the "time-space" components of left hand side Eq. ([69]) . 



We can thus solve Eq. (|119[) for hoi, 



hoi — 



8nG{0\ o - /3| ) 



/9i 



4(1 - 87tG/?i ) + 167rG£ - B a C a 



Vi 



16ttGE - B C - - ) H ,; 



(121) 



where we have used that the superpotential \ satisfies 
X,oi = Vi-Wi. 



5.3 Solving g 00 to 0(2) 

A full analysis of the PPN parameters requires that we 
solve for the "time-time" component of Eq. (f6"9")> to 0(2). 
To 0(2) we have 



2U + 2U Z - 8<7> 2 - 2B CoX,oo (122) 



Roi = —^Oijj + ^(1 + 2BC ) 
and the "time-space" components of T^, 



T, 



(m) 



0; 



-o v- T [<k) - 



and 



rp(u) _ j 7 

± 0i — ^0i,0 > 

where Joi,o is given by Eq. (|107[) and 



(114) 

(115) 
(116) 



Ji 



-(010 + P20)u 3 tij ~ PmUjj - (010 - ^3o)/lo[,',i] 



'2(^10 + 03a)hij,Oj — 2^ 20 ^ii> 0i • 



(117) 



Using Eqs. (I112[) and (|113p . and our previous results, we 
obtain 



T («) (lio-fffoh 



20 



10 



10 



(/3i 2 o-/3 3 2 o)(3-2BoG ) 



4/3 



10 



X,oijj -(118) 



Here, we have assumed that B is a constant parameter, 
B — > i?o, which will be determined below. 

By solving the "time-space" components of the field 
equation, we obtain 



8ttG(/4 - g|p) 
0io 



V 



l6nGE - B Co ~ - ) Xfiijj , (H9) 



where 



(/3 2 o-/?3 2 o)(3-2i3oCo) 
4/3io 



(120) 



where we have defined h o = goo + 1 — 2U. To 0(2), we 
also have 



r (m) 



J °° - 2 



p m (l + 77 + - 217), T^ m) = p m UjUj + 3p ri 



j,(l)\2 T W 



/,( 1 )^2 



) 2 , 



W =(3i(V uy + 2u ViJ \ 
=/3 1 (V ^) 2 -2V J ll . 



(123) 

Using the covariant derivative Eqs. (|89|) and (|90p . the com- 
ponents of energy-momentum tensor for the vector field 
become 



T, 



(u) 



and 



00 



-01 



2U +-U 2 -9$ 2 
2 



-0 1O [3-2C O (1+B O )] X S)0ii , 



(124) 



■(fiia + 3/3 20 + /3 30 )(3 - 2C ) Xmi ■ (125) 



From Eq. (|123[) one can evaluate the "time-time" compo- 
nents of the right hand side of Eq. 



T, 



,(m) 1 rnkklf.i/f/lMifnMl 



1 - 8irGf3 



10 



BttG 



(U + 2<7>i - 2<2> 2 + ®3 + 3^4) ,(126) 



- \goo9,uT^ = \(T$ + T^) = 0, (127) 



and 



T oo - -j^oog^v 



1 



rp(u)nu _ 1 (rp[u) . rp(u)\ 



/ l0 o + 2U + 2U Z - 84-2 
-{20x0 + 302 + /3so)(3 - 2G )]x,oori. 



-[275 G 



(128) 
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Using the "time-time" component of Eq. (|69|) to 0(2) and 
combining Eqs. (fT22]) . (fT26]l . (fT27]) and (fT28|) . we can ob- 
tain 

h 0Q = -2U 2 + 4<Z>i + 4<Z> 2 + 2<2> 3 + 6<2> 4 

+QoX,oo , (129) 



where 



Qo = 



16ttGC q 
1 - 8ttG/3 



10 



/3io + 2/? 3 o + 



8ttG/3 



10 



8ttG 



B ) (.130) 



The last term of Eq. (|129p is a new PPN parameter. How- 
ever, there is no need to introduce any additional PPN 
parameters when we move into the standard gauge by 
choosing Bq such that Qq = 0. We have 



B 



8irG(f3 w + 2/3 30 ) 



1 - 8nG(3 



10 



(131) 



We have all the necessary ingredients to read off the 
PPN parameters. Let us begin by writing the full metric, 
we have 

5oo = -1 + 2U - 2U 2 + 4<Z> 2 

+4<?i + 2<P 3 + 60 4 , 
g i:j = (1 + 2U)5 l0 , 

#10 



.9o i 



4(1 - 8itG[3 10 ) + 16itGE - B„C - - 



(l6irGE - B C - i J Wi 



Vi 

(132) 



By comparing Eqs. ([74]) and (|132|) . the PPN parameters 
are given by 



Ippn = 1, Pppn = 1, X — Ci = C2 = C3 = C4 = a 3 = 0, 
/ 1 



^30 



«2 = 



/? 10 (l-8^G/3 10 ) + 8^G/3 2 ' 
(0 W + 2/3 30 ) 2 
V8^G; z (1-8^G/3 10 )(/3io+/52o + /3 3 o) 
6/3io/9 3 o + 3^0 [1 - 8ttG(I3 10 + 2/j 30 )] 
(1 - 8ttG(3 w ) [/3io(1 - 8ttG/3 10 ) + 8ttG/3 3 2 ] 

4/3 3 2 [l-27rG(/3 10 -6/3 30 )] 
(1 - 8vrG/3 10 ) [/3io(l - 8ttG/3 10 ) + 8ttG/3 3 2 ] " 



(133) 



Notice that the values of "/ppn and (3ppn are the same 
as in GR. jppn = 1 implies an identical predicted deflec- 
tion of light about the Sun as well as identical predic- 
tions for radar echo delay. The parameter (3ppn enters 
into the expression for anomalous relativistic precession 
of planetary orbits, but the strongest experimental limit 
is provided by the lunar laser ranging test. The other non- 
vanishing ones are ot\ and ct2 which are the effect of the 
preferred frame. 



6 Conclusions 

In this paper, we have investigated the cosmological evo- 
lution of an interacting scalar field model in which the 
scalar field has an interaction with the background mat- 
ter via Lorentz violation. We propose a model of interac- 
tion, specifically Q m = —f3p m /(3 in which the interaction 
is mediated by the slope of the effective coupling (3. The 
equation of state parameter of the scalar field is expressed 
by eq. ([34|) as a candidate of dark energy. The important 
role of the model is played by the effective coupling in 
the transition era from the matter dominated to scalar 
field dominated, which leads to an accelerating universe. 
The model also predicts a constant fraction of dark en- 
ergy to dark matter in the future and hence solves the 
coincidence problem. This is a profitable support to the 
effective coupling. As a cosmological implication, the dy- 
namic of the effective gravitational constant is determined 
by the effective coupling and allows one to test the Lorentz 
violating scalar-vector-tensor theory of gravity using the 
Snla data. We have studied how a varying G or a effective 
coupling could modify the evolution of the Hubble param- 
eter which is deviated by the term of j3~ 2 . For a simple 
polynomial (3{z) = (3q(1 + £z 2 ) ansatz, the best fit values 
are xL™ = 195.68, ( = -0.33, and O m0 = 0.24 for the 
modified /LCDM model and xL™ = 195.71, C = -0.29, 
and ujfj, — —1.13 for the modified quintessence model. 

We also have presented the IPPN parameters of the 
theory. Our result strongly depends on the two assump- 
tions: the scalar field satisfies the extremum condition, 
and the coupling functions (3\ and P2 are constant. The 
first assumption stems from the fact that the potential 
will play the role of an effective cosmological constant if 
the theory is to account for the late time accelerated ex- 
pansion of the universe. The second one is to obtain the 
real velocity of the gravitational waves. Up to IPPN ap- 
proximation, important to note that it is impossible to ob- 
tain Gn — G cosmo , although fix and 02 are constant. Here 
Gcosmo = G (e) is given by Eq. ((50j) . G cosmo is still dynamic 
in this case because of P2, while Gn is always constant 
(|99|) . In the post-Newtonian approximation, the effect of 
non-constant coupling function 02 appears in 0(2.5) or 
higher. 

Of course, there are many remaining works to make 
this scenario more concrete which is beyond the main aim 
of the present work. For instance, non-linear coupling or 
more complicate functions are also possible. In the present 
work we have considered to the case of instantaneous crit- 
ical points, where Ai and A2 are the constant parameters. 
However, in order to construct viable model, Ai and A2 
should satisfy the condition Eq. (I37|) . We need to solve the 
dynamical system Eq. (|36[) . Then, Ai and A2 are dynam- 
ically changing quantities. For example, the critical point 
D in Table Q] becomes x(N) = y / 372/(Ai(iV) + X 2 (N)), 
y(N) = x(N), z(N) = 0, and u(N) = where N = Ina. 
As a consequence we obtain the running critical point ac- 
cording to the changing of Xi(N) and X2(N). 

Another important aspect is that the coupling interac- 
tion between the scalar field and matter fields affects only 
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the solution for the scalar field. <jP^ represents the local 
deviation from <pQ, which vanishes far from the local sys- 
tem. In conclusion, one can reasonably state that the our 
gravity model could be a viable candidate theory, even in 
the PPN approximation. It cannot be a priori excluded at 
Solar system scales. 
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